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Abstract 

Background: Previous studies have shown that the activity of the amygdala is elevated in people experiencing clinical and 
subclinical levels of anxiety and depression (negative affect). It has been proposed that a reduction in inhibitory input to the 
amygdala from the prefrontal cortex and resultant over-activity of the amygdala underlies this association. Prior studies 
have found relationships between negative affect and 1) amygdala over-activity and 2) reduced amygdala-prefrontal 
connectivity. However, it is not known whether elevated amygdala activity is associated with decreased amygdala- 
prefrontal connectivity during negative affect states. 

Methods: Here we used resting-state arterial spin labeling (ASL) and blood oxygenation level dependent (BOLD) functional 
magnetic resonance imaging (fMRI) in combination to test this model, measuring the activity (regional cerebral blood flow, 
rCBF) and functional connectivity (correlated fluctuations in the BOLD signal) of one subregion of the amygdala with strong 
connections with the prefrontal cortex, the basolateral nucleus (BLA), and subsyndromal anxiety levels in 38 healthy 
subjects. 

Results: BLA rCBF was strongly correlated with anxiety levels. Moreover, both BLA rCBF and anxiety were inversely 
correlated with the strength of the functional coupling of the BLA with the caudal ventromedial prefrontal cortex. Lastly, 
BLA perfusion was found to be a mediator of the relationship between BLA-prefrontal connectivity and anxiety. 

Conc/usions:J\r\ese results show that both perfusion of the BLA and a measure of its functional coupling with the prefrontal 
cortex directly index anxiety levels in healthy subjects, and that low BLA-prefrontal connectivity may lead to increased BLA 
activity and resulting anxiety. Thus, these data provide key evidence for an often-cited circuitry model of negative affect, 
using a novel, multi-modal imaging approach. 
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Introduction 

One commonly cited model of mental illness is derived from the 
writings of the neurologist Hughlings Jackson, who observed that 
frontal lobe lesions were associated with 'release behaviors', i.e., 
disinhibition of suppressed, primitive reflexes that are generated 
subcorticaUy [1,2]. Based on Jackson's ideas, theories about 
circuitry abnormalities underlying a diverse set of disorders, 
including depression [3,4], Post-Traumatic Stress Disorder [5], 
and schizophrenia [6,7], have suggested that disruptions of 
projections from the prefrontal cortex to subcortical structures 
lead to a loss of inhibitory control and over-activity of these 
structures, with the subsequent emergence of symptoms. Neuro- 
imaging research over the past several decades has generated 
support for this model [8-1 1], most consistently with demonstra- 



tions of associations between increases in negative affect (i.e., 
anxiety, depression, or subsyndromal expressions of those symp- 
toms [12]) and 1) increased responses [13] or resting activity 
[14,15] of the amygdala, or 2) decreased strength of the 
connections between the prefrontal cortex and amygdala [IB- 
IS]. However, it has remained unclear whether decreased 
prefrontal-amygdala connectivity is associated with amygdala 
over-activity during these states. 

A novel approach to testing this model is to combine two MRI- 
based methods that measure "resting-state" brain activity: 1) 
arterial spin labeling (ASL) functional MRI (fMRI), which uses 
magnetically-labeled arterial blood water as an endogenous tracer, 
providing a quantitative measure of regional cerebral blood flow 
(rCBF), and 2) blood oxygenation level-dependent (BOLD) 
functional connectivity fMRI, which measures correlations in 
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spontaneous, low frequency fluctuations in tlie BOLD signal 
between diflerent brain regions [19]. ASL fMRI has been 
validated by comparisons with positron emission tomography 
measurements of rCBF at rest [20]. It is particularly well-suited for 
the quantitative measurement of tonic neural processes, since it is 
relatively free of the slow drift associated with BOLD contrast 
imaging, and has good reliability [20,21]. Resting-state BOLD 
functional connectivity, on the other hand, is now widely used as a 
measure of the degree of functional coupling among distinct brain 
regions comprising a network [22]. 

Previous studies using resting-state BOLD data to measure 
functional connectivity have shown that the basolateral nucleus of 
the amygdala (BLA) in humans displays a high degree of 
functional correlation at rest with the ventral medial prefrontal 
cortex (mPFC), whereas the central nucleus of the human 
amygdala shows little connectivity with the mPFC [23,24], similar 
to the pattern of amygdala projections found in non-human 
primates [25,26]. In the present study, we first tested whether ASL 
fMRI could detect a relationship between one type of negative 
affect, current anxiety levels, and BLA activity. We then used ASL 
and BOLD fMRI in combination to test the hypothesis that 
elevated BLA rCBF is associated vnth increased anxiety and 
reduced BLA-prefrontal connectivity. 

Materials and Methods 

Ethics Statement 

Written informed consent was obtained from all subjects prior 
to enrollment, in accordance with the guidelines of the Partners 
Healthcare Institutional Review Board, which approved this 
study. 

Subjects and Anxiety Measure 

38 healthy adults ( 1 1 females; see Table 1 for demographic 
information) were recruited via advertisement. Subjects were 
excluded if they had a current psychiatric disorder or history of 
one, including a history of substance abuse or dependence, or who 
had used illicit substances within the three months prior to the 
study, as determined by the Structured Clinical Interview for 
DSM-IV (SCID) [27]. Subjects were also excluded if they had ever 
used psychotropic medications, had any contraindications for 
MRI scanning (metal implants, claustrophobia), or had a 
neurologic disorder or a history of head trauma with neurologic 
sequelae. Four subjects did not have resting BOLD scans and thus 
were not included in the functional connectivity MRI analyses. 
Four subjects did not have ASL scans and thus were not included 
in the ASL analyses. One additional subject was excluded from the 
ASL analyses due to poor data quality (movement >3 mm). The 
Spielberger State and Trait Anxiety Inventory (STAI) [28] was 
administered on the same day that the scanning occurred. Factor 
analyses have shown that the STAI measures symptoms typically 
associated with depression (e.g., anhedonia), as well as those 
associated with anxiety; thus it appears to index the broad 
construct of negative affect, rather than anxiety only [29—34]. 
Consistent with this, it has been found that scores on the STAI are 
equally or more highly correlated with measures of depression 
than with scores on other anxiety scales [32]. However, since we 
used just one measure of negative affect (the STAI state measure), 
we have referred to the data collected with the STAI here as 
"anxiety levels". 

Image Acquisition 

Imaging took place in a 3-Tesla Siemens TIM Trio MRI 
scanner (Siemens Medical Systems, Erlangen, Germany) using a 



12-channel head coU. Subjects underwent two conventional high- 
resolution three-dimensional structural Tl magnetization pre- 
pared rapid acquisition gradient-echo (MPRAGE) scans (8 min 
7 sec, 128 sagittal slices, 1.33-mm thickness, repetition time 
(TR) = 2.530 ms, echo time (TE) = 3.39 ms, flip angle = 7°, reso- 
lution 1.3x1x1.3 mm). 

One 6-min-2-sec pulsed ASL scan was collected using the 
"PICORE-Q2TIPS" MRI labeling method [35], with the 
following parameters: TR — 3000 ms, TE =13 ms, TI 1 = 700 ms, 
TI2= 1700 ms, voxel size = 3.515 x3. 515 x6. 25 mm, number of 
slices = 16 (5 mm thickness with a 1.25 mm gap), FOV = 22.5 cm, 
matrix = 64x64, flip angle = 90 degrees. 'Tag' images were 
acquired by labeling a thick inversion slab (110 mm), proximal 
to the imaging slices (gap = 2 1.2 mm). 57 'Control' images were 
acquired interleaved with 57 tag images by applying an off- 
resonance inversion pulse without any spatial encoding gradient. 
At the beginning of each ASL scan, an Mq scan (i.e., the 
longitudinal magnetization of fully relaxed tissue) was acquired for 
resting regional cerebral blood flow (rCBF) quantification purposes 
[36]. A delay time of 1 second was inserted between the end of 
labeling pulses and image acquisition to reduce transit-related 
effects. 

One 6-min-12-sec resting BOLD scan (TR= 3000 ms; 
TE = 30 ms; flip angle = 85 degrees; 47x3 mm thick slices, 
3x3x3 mm in-plane resolution) was also acquired. 

During both the ASL and BOLD scans, subjects were instructed 
to keep their eyes open, bfink normally, and look at a blank screen. 

Regions-Of-lnterest 

Regional CBF levels were extracted from the basolateral 
nucleus of the amygdala (BLA) using an anatomically-defined 
BLA region-of-interest (ROI). The BLA ROIs were derived from 
the Jeulich Histological Adas [37], as in previous studies of 
amygdala functional connectivity [23,24]. The ROIs were 
comprised of voxels with £50% probability of being located in 
the BLA. Voxels near a subregion border were included if they 
were more likely to be located in the BLA than in other amygdala 
subregions [23]. For the ASL analyses, these ROIs were then 
automatically mapped into each individual subject's native 
anatomical scan. To maximize power, a seed combining the left 
and right BLA ROIs was created and used as the seed for the 
functional connectivity fMRI analyses. 

ASL fMRI Analysis 

ASL data analyses were performed using a combination of FSL 
(FMRIB's Software Library, www.fmrib.ox.ac.uk/fsl) [38] and 
FreeSurfer (http://surfer.nmr.mgh.harvard.edu/) tools. Tlu- 'tag', 
'control', and Mo scans were first motion-corrected using 
MCFLIRT [39] and r(-gist(;red to the high resolution anatomical 
images using FreeSurfer's bbregsiter tool [40]. Then tag and 
control scans were surround subtracted (i.e. given each tagx, 
[(controlx-i+controlx+i)/2 — tagx]) to generate perfusion-weight- 
ed images. All of the perfusion-weighted maps were then averaged 
and scaled by a factor proportional to the Mg scan to obtain rCBF 
maps in absolute values (ml/100 g of tissue/min) [36]. Group-level 
analyses were conducted in the volume, and group rCBF maps 
were smoothed with a 3D isotroj)ic, 7.02 mm full-width half- 
maximum (FWHM) Gaussian kernel (set to 2*voxel size). Regional 
CBF values were extracted from the BLA ROIs mapped into each 
subject's individual anatomical scan. Additionally, whole-brain 
voxel-wise analyses using STAI state anxiety score as a regressor 
were conducted, in order to determine whether there were 
significant correlations between anxiety and rCBF of regions other 
than the amygdala. In this analysis, the Gaussian Random Fields 
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Table 1. Participant demographics. 





Mean (SD) 


Range (Min, Max) 


Age (years) 


29.7 (9.9) 


19, 53 


Education Level {years completed} 


15.7 (1.5) 


12, 19 


Verbal IQ^ 


111.2 (6.1) 


97.6, 125.1 


Anxiety (state)'' 


25.2 (4.4) 


20, 37 


Anxiety (trait)*" 


27.4 (7.1) 


20, 47 


Socioeconomic Status*^ 


47.3 (8.6) 


32, 66 


Depression'^ 


1.1 (1.7) 


0, 6 


Caucasian/Non-Hlspanic (%) 


78.9 




Caucasian/Hispanic (%} 


7.9 




Black/African-American (%) 


7.9 




Asian (%) 


5.3 





Participant demographics are listed. Values given are Mean (Standard Deviation) unless otherwise specified. 

^Verbal IQ was measured using the North American Adult Reading Test [74]. 

""Anxiety was measured using the Spielberger State and Trait Anxiety Inventory [28]. 

'^Socioeconomic Status was measured using the Holllngshead Scale [75]. 

'^Depression was measured using the Beck Depression Inventory [76]. 

doi:1 0.1 371 /journa!.pone.0097466.t001 



(GRF) method was used to correct for multiple comparisons [41], 
with a voxel-level height threshold of p<.0005. Clusters that 
showed cluster-wise significance (cluster-level threshold of p<.05, 
GRF whole brain corrected) are reported in the text and in 
Table 2. 

Functional Connectivity fMRI Analysis 

Standardized preprocessing techniques [42] were used to 
selectively capture variance in the BOLD signal corresponding 
to low-frequency (<.08 Hz) fluctuations in neural activity. Data 
were smoothed using a 7 mm FWHM Gaussian kernel. Nuisance 
regressors, including the six parameters computed from the rigid- 
body motion correction, the averaged signal within a ventricular 
ROI, a region within the deep white matter, and the signal 



averaged over the whole brain, were used to remove systematic 
variance associated with these variables. The first temporal 
derivative of each regressor was also included to account for 
temporal shifts in the BOLD signal. Measures of slice-based and 
voxel-based signal-to-noise ratios (SNR) and the following motion 
parameters [43] were calculated for each BOLD scan to test if 
these parameters correlated with anxiety or perfusion levels: 1) 
Mean relative translations in 3D space (each volume compared to 
the previous timepoint), 2) Mean absolute translations in 3D space 
(each volume compared to the first timepoint), 3) Number of 
micro-movements (the number of relative translations in 3D space 
>.5 mm), and 4) Mean absolute rotations in the x, y, and z 
directions. 



Table 2. Regions showing a correlation between perfusion and anxiety levels. 



Region 


BA 


HemI 


Tal (x, y, z) 


Size (mm^) 


Z 


Superior frontal gyri 


9 


R 


9, 50, 25 


7811 


5.0 


Superior parietal gyrus, precuneus 


7 


R 


23, -72, 34 


1197 


4.5 




7 


R 


14, -52, 50 


1771 


4.4 


Thalamus, hypothalamus, posterior hippocampus 




L 


-11, -20, -1 


6714 


4.5 






R 


13, -21, -1 


3002 


4.2 


Putamen (posterior), hippocampus, amygdala 




L 


-34, -17,-7 


3283 


4.5 


Putamen (anterior) 




L 


-11, 19, -11 


1769 


4.5 


Postcentral gyrus 


41 


L 


-48, -19, 23 


1145 


4.4 




43 


R 


60, -11, 13 


948 


4.0 


Middle frontal gyrus 


46 


R 


30, 45, 24 


2286 


4.3 


Posterior cingulate gyrus 


23 


L 


-5, -43, 24 


3504 


4.3 


Rostral anterior cingulate, medial frontal gyri 


32/9 


L 


-7, 47, 11 


525 


4.1 


Hippocampus, putamen 




R 


41, -19, -9 


2094 


4.0 



Results of a whole-brain voxel-wise regression analysis of the regional cerebral blood flow (rCBF) data, using anxiety levels as a regressor, are listed. Sites which showed 
a significant positive correlation between rCBF and anxiety levels are listed below. Also, see Figure 2. BA = Brodmann Area; Hemi=; hemisphere; Tal = Talaraich 
coordinates. 

doirl 0.1 371/journal.pone.0097466.t002 
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To create whole-brain correlation images, the averaged time 
series across all voxels comprising the seed (the bilateral BLA ROI) 
was used as the variable of interest in a linear regression with the 
time series corresponding to each voxel across the brain. All 
statistical analyses of correlation data were performed on Fisher z 
transforms [44]. Only clusters that reached whole-brain family- 
wise error (FWE) significance with a height threshold of p<.05 
and size threshold of S5 voxels are reported (Table 2). 

Whole-brain random effects multiple regression analyses using 
STAI state anxiety score or rCBF of the BLA as effects of interest 
were conducted to identify locations showing correlations between 
the strength of their connectivity with the seed and the anxiety or 
BLA rCBF measures. Analyses were restricted to regions of the 
brain that showed significant connectivity with the seed region 
(either positive or negati\ K). Based on previous work showing that 
the BLA is highly anatomically and functionally connected with 
the ventral mPFC (see Introduction), and a recent meta-analytic 
factor analysis suggesting that the ventrocaudal mPFC is primarily 
involved in 'affect generation' [45], we chose the ventral mPFC as 
the a priori search area for the regression analyses. The search area 
was confined to a ventral mPFC mask (39,192 mm^) consisting of 
voxels showing significant BLA-mPFC functional coupling within 
arbitrary boundaries containing the ventral mPFC region: x- 
coordinate £ | 10 1 , y-coordinate ^0, and z-coordinate £10. Thus 
this independent, a priori mask was used in order to restrict the 
subsequent analyses to voxels within the ventral mPFC showing 
functional connectivity with the BLA, prior to conducting the 
correlational analyses of interest. Sites showing correlations 
between BLA-mPFC connectivity and the regression \'ariable 
(anxiety or rCBF) were deemed significant if they met a cluster- 
wise correction (FWE, p<.05) within the ventral mPFC. All 
second-level analyses were conducted using SPM5 (http://www. 
fil.ion.ucl.ac.uk/ spm). 

The location of significant clusters was determined using the 
Talairach and Tournoux Stereotaxic Atlas [46] and confirmed 
with the Wake Forrest University (WFU) PickAtlas [47]. 

Secondary Mediation Analyses 

A mediation analysis was then conducted to assess the 
relationships among the results of our regression analyses (i.e., 
the direction of those effects, with two possible outcomes or models 
of the relationships, see below). The relationship between BLA- 
mPFC connectivity strength, BLA perfusion levels, and anxiety 
scores was assessed using the Indirect Mediation Analysis script 
created by Preacher and Hayes for SPSS [48]. Non-parametric 
bias-corrected bootstrapping (with 10,000 bootstrap samples) was 
used to derive 95% Confidence Intervals of the indirect effects. 

Results 

Anxiety 

The mean (± standard deviation) state anxiety score of this 
sample (25.2 ±4.4; median = 24, range 20-37, with a possible 
range of 20-80) was well below means typically reported in 
patients with anxiety and other psychiatric disorders [28,49]; 
scores on the STAI-state of greater than 40 are considered 
clinically significant [49]. Thus, all subjects had STAI-state scores 
in the healthy to subsyndromal range. See Table 1 for additional 
details about the subjects. 

ASL Analysis 

Perfusion of both the right and left BLA was significandy 
correlated with anxiety levels (both right and left BLA: r = .54, p< 
.002) (Figure lA and IB). This finding was then confirmed in a 



follow-up, voxel-wise, whole brain-corrected regression analysis, 
which revealed strong correlations between anxiety and perfusion 
of voxels within the right (x, y and z Talairach coordinates of peak 
voxel: 33, -5, -13; p<5xl0"*) and left (-22, -10, -22; p< 
1x10 ^) amygdala, the majority of which fell within the 
boundaries of the BLA ROIs (Figure IC). In addition to the 
amygdala, the whole brain regression revealed that anxiety was 
also significandy correlated with rCBF of the superior frontal gyri, 
posterior cingulate cortex (PCC), thalamus, and putamen 
(Figure 2A and 2B), as well as the superior parietal gyrus and 
precuneus, postcentral gyrus, middle frontal gyrus, rostral anterior 
cingulate and medial frontal gyri, and hippocampus (Table 2). 

Resting-state Functional Connectivity Analysis 

An average map of the functional connectivity of the BLA 
confirmed [23,24] that BLA activity is significantly coupled with 
the activity of the ventral mPFC, hypothalamus and brainstem 
(Figure 3A; Table 3). 

Whole Brain Regression Analyses 

Whole brain regression analyses revealed that anxiety was 
significandy inversely correlated with connecti\'ity between the 
BLA and several regions of the mPFC, including the subgcnual 
cingulate gyrus (SGC) (extending to the adjacent hypothalamus; 4, 
2, —7; Pcorrected <-02; 3392 mm^; Z = 4.4), orbitofrontal cortex, 
and superior frontal g^Tus (Figure 3B). Moreover, BLA rCBF was 
also inversely correlated with connectivity between the BLA and 
the mPFC, including the SGC (extending to the adjacent 
hypothalamus; 2, 4, -4; Pcorrected <-04; 3440 mm^; Z = 4.2) and 
medial frontal cortex (Figure 3D). These peaks were located in the 
posterior-most portion of the SGC bordering on the hypothalamus 
[46,47]. In summary, lower BLA-mPFC functional connectivity 
was associated with both greater BLA rCBF and greater anxiety. 

None of the motion parameters (relative or absolute head 
displacement, head rotations, or micro-movements) or SNR 
measures correlated with anxiety or perfusion (all rs<.32, all 
ps>.10). 

Mediation Analysis 

A follow-up, secondary mediation analysis was then conducted 
to assess the relationship among the findings of the regression 
analyses conducted above, to distinguish between the following 
two models (also see Figure 4): 1) BLA perfusion mediates the 
relationship between BLA-mPFC connectivity and anxiety (Model 
7^1), or 2) BLA-mPFC connectivity mediates the relationship 
between BLA perfusion and anxiety (Model #2). 

To conduct this analysis, BLA-mPFC connectivity levels that 
had been found to correlate with anxiety or BLA perfusion were 
extracted for each subject in the following manner: The time 
course of mPFC resting-state BOLD activity was extracted from 
3 mm-radii spheres centered on the locations of the peak 
correlations with anxiety levels (for the test of Model ^\) or 
BLA perfusion (for the test of Model #2), and then correlated with 
the time-course extracted from the anatomical BLA seed. These 
correlations were then transformed into z-scores for each subject, 
reflecting BLA-mPFC connectivity strength. 

Testing of Model jj^l (Figure 4A) revealed a significant direct 
effect of BLA-mPFC connectivity on perfusion levels (path a: — 
21.88±5.48 [unstandardized coefficient ± standard error], t= — 
3.99, p<.001), a significant direct effect of perfusion levels on 
anxiety (path b: .23±.09, t = 2.46, p<.05), and a significant total 
effect of BLiA-mPFC connectivity on anxiety (path c: — 
11.06±2.86, t=— 3.87, p<.001). However, when perfusion was 
taken into account as a mediating variable, the direct effect of 
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Anxiety Anxiety p < .005 



Figure 1. Perfusion of the basoiateral nucleus of tKie amygdala (BLA) is correlated with anxiety levels. Significant correlations were 
found between anxiety levels and perfusion of the left (A) and right (B) basoiateral amygdala (BLA), as defined using anatomical regions-of-interest. 
These findings were then confirmed in a voxel-wise, whole brain regression analysis (C). In C, the BLA regions-of-interest are outlined in blue; the 
voxel-level display threshold is p<.005 (showing only clusters surviving whole-brain correction, see Methods). Clusters that showed cluster-wise 
significance (p<.05, whole brain corrected) are reported in the text and in Table 2. R, right. 
doi:10.1371/journal.pone.0097466.g001 



BLA-mPFC connectivity on anxiety was no longer significant 
(path c': — 6.10±3.31, t= — 1.84, p = .08), indicating a specific 
mediated effect. Additionally, the bias-corrected bootstrapping test 
of the indirect effects of BLA-mPFC connectivity on anxiety 
through BLA perfusion as a mediator (path ab: — 5.44±3.05) 
reached significance, with the 95% confidence intervals (—13.38, 
— 1.27) indicating a specific mediated effect. In other words, BLA 
perfusion levels accounted for a significant portion of the 
relationship between BLA-mPFC connectivity and anxiety. 

Testing of Model #2 (Figure 4B) revealed a significant direct 
effect of perfusion levels on connectivity (path a: — .02 ± .003, t = — 
6.04, p<l xlO a significant total effect of perfusion levels on 
anxiety (path c: .33±.08, t=4.33, p<.OOI), but a non-significant 
direct effect of connectivity on anxiety (path b: 1.79±4.58, t= .37, 
p = .71). Additionally, the direct effect of perfusion levels on 
anxiety remained significant after taking into account BLA-mPFC 
connectivity strength (path c': .36±.12, t = 3.06, p<.01), indicat- 
ing that the relationship between BLA perfusion and anxiety is not 
mediated by BLA-mPFC connectivity. Furthermore, the bias- 
corrected bootstrapping test of the indirect effects of BLA 
perfusion on anxiety through BLA-mPFC connectivity as a 



mediator (path ab: —.03 ±.10) did not reach significance, with 
the 95% confidence intervals (—.21, .18) failing to indicate a 
specific mediated effect. In other words, BLA-mPFC connectivity 
failed to account for a statistically significant portion of the 
relationship between BLA perfusion levels and anxiety. In sum, 
these mediation results support Model #1, that reduced BLA- 
mPFC connectivity leads to elevated BLA perfusion and subse- 
quent increases in anxiety, rather than the alternative possibility 
that elevated BLA perfusion leads to poor BLA-mPFC connectiv- 
ity and increases in anxiety (Model 41^2) . 

Discussion 

The goal of this study was two-fold. First, we aimed to 
determine whether ASL fMRI could detect changes in amygdala 
activity related to variation in negative affect in healthy subjects. 
We found strong correlations between the perfusion of the BLA, as 
well as the hippocampus, putamen, and multiple midline frontal 
and parietal cortical areas, and anxiety levels, suggesting that ASL 
is a sensitive tool for measuring this relationship. Second, we 
sought to test a widely-cited model of psychopathology-the 
proposal that certain symptoms, including clinical anxiety and 




Figure 2. Perfusion of a distributed networic of regions outside of the amygdala is also correlated with anxiety levels. A voxel-wise 
whole brain regression analysis revealed that, in addition to the basoiateral amygdala (BLA), perfusion of the superior frontal gyri and posterior 
cingulate cortex (A), and anterior putamen (B), among other regions (see Table 2), were significantly correlated with anxiety levels. Whole-brain 
corrected results (see Methods) are displayed here using a voxel-level threshold of p<.005. Clusters that showed cluster-wise significance (p<.05, 
whole brain corrected) are reported in the text and in Table 2. R, right; PCC, posterior cingulate cortex; SFG, superior frontal gyri. 
doi:1 0.1 371/journal.pone.0097466.g002 
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Table 3. Basolateral amygdala functional connectivity. 



Region 


BA 


Heml 


Tal (x, y, z) 


Size (mm') 


Z 


Medial orbitofrontal gyrus 


11 


L 


-6, 19, -18 


2160 


6.4 




n 


L 


-6, 58, -11 


40 


5.4 


Anterior medial orbitofrontal gyrus, frontal pole 


10 




0, 66, 3 


152 


5.7 


Hippocampus 




R 


20, -33, 0 


72 


5.6 


Insula 




R 


32, -7, 10 


40 


5.5 


Lateral orbitofrontal gyrus 


n 


R 


28, 28, - 1 7 


56 


5.3 


Inferior parietal gyrus 


7 


R 


50, -45, 35 


1816 


6.5 


Precuneus, posterior cingulate gyrus 


7 


R 


4, -69, 51 


4488 


6.3 




7/31 


R 


4, -38, 46 


264 


5.5 




23 


R 


4, -28, 25 


192 


5.5 


Inferior middle frontal gyri 
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Areas of the brain showing significant functional connectivity with the basolateral amygdala {BLA} are listed. Clusters that are unshaded are those with positive 
functional coupling with the BLA, whereas clusters that are shaded grey are those showing negative functional coupling (inverse or anti-correlations) with the BLA 
(following global mean regression). Sites of connectivity within or abutting the BLA are not listed because of the difficulty of interpreting these findings. Also see 
Figure 3. BA^Brodmann Area; Hemi = hemisphere; Tal =Talaraich coordinates. 
doi:l 0.1 371 /journal.pone.0097466.t003 



depression (which are tliought to exist on a continuum with 
subsyndromal, everyday anxiety and other negatively-valenced 
affective states [50-52]), arise from dysregulation and resulting 
over-activity of subcortical structures such as the amygdala. This 
over-activity has been proposed to arise from a reduction in 
inhibitory input from the prefrontal cortex. Here, using advanced 
MRI-based techniques and data collected in healthy subjects who 
have not been exposed to potentially confounding treatments and 
chronic illness, we detected significant associations among all three 
components of this model, reduced prefrontal-subcortical (BLA) 
connectivity, subcortical (BLA) over-activity, and anxiety, in the 
same subjects. 

Our finding of a strong correlation between anxiety levels and 
BLA perfusion are consistent with prior Positron Emission 
Tomography findings of associations between depression severity 
and the magnitude of amygdala blood flow [14] or metabolism 
[15] in patients with Major Depressive Disorder, and associations 
between levels of negative affect and amygdala blood flow in 
healthy subjects [53]. Also, task-based fMRI studies have found 
that the magnitude of amygdala BOLD responses is correlated 
with anxiety levels in healthy subjects [54-56]. The present 
fmdings reveal that ASL fMRI can also be used to detect and 
measure this relationship, providing a means of measuring the 
association that does not require exposure to radiotracers or a 
capacity to attend and respond to experimental stimuli. 

Secondly, the results of the whole-brain regression analyses 
conducted here replicate and extend previous work showing that 
subsyndromal negative affect (measured with the same scale used 
here or with measures of overlapping constructs such as 
neuroticism) is predicted by the strength of amygdala-mPFC 



connectivity. Fezawas and colleagues first reported that amygdala- 
mPFC connectivity, as measured by task-elicited BOLD responses, 
is inversely correlated with temperamental negative affect [16]. 
This finding was later rephcated using resting-state connectivity 
methods [17]. Also, Kim and Whalen used diffusion tensor 
imaging to show that anxiety levels are inversely correlated with 
the structural integrity of a white matter pathway overlapping with 
the expected trajectory of the uncinate fasciculus [18], a fiber tract 
containing projections between the amygdala and the subgenual 
cingulate and orbitofrontal gyri [57,58]. This association may have 
other structural manifestations; for instance, negative affect levels 
have been correlated with increased amygdala volumes and 
decreased thickness of the left subgenual and rostral anterior 
cingulate cortices in healthy young adults [59]. Thus, our finding 
of an association between anxiety and amygdala-mFFC connec- 
tivity extends this body of previous work. 

The present study shows for the first time that amygdala hyper- 
perfusion is linked to both anxiety and poor connectivity between 
the BLA and the caudal-ventral mFFC (BA25/subgenual cingu- 
late cortex (SGC)) in the same subjects. Because conclusions 
cannot be drawn about the direction of these associations using 
regression alone, we used a mediation approach to determine 
whether, in our subjects, reduced connectivity strength within the 
BLA-mPFC pathway could be hnked causally to elevated BLA 
perfusion and resulting anxiety. These analyses showed that the 
inverse relationship between BLA-mPFC connectivity and anxiety 
was significantly mediated by BLA perfusion levels, i.e., that 
decreased connectivity between the BLA and ventral mPFG 
contributes to BLA overactivity, resulting in higher anxiety levels. 
The alternative model, that BLA-mPFC connectivity mediates the 
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A. BLA Functional Connectivity (FC) 





Anxiety 

D. BLA FC - BLA rCBF Correlations E. 




Figure 3. Functional connectivity between the BLA and mPFC is inversely correlated with BLA perfusion and anxiety levels. An 

average map of basolateral amygdala (BLA) functional connectivity is shown in A. A whole-brain voxel-wise regression revealed that the strength of 
connectivity between the BLA and mPFC was negatively correlated with both: anxiety levels (B, C) and BLA perfusion (D, E). In A, B, and D, voxels with 
positive connectivity with the BLA (A) or showing positive correlations between their connectivity with the BLA and anxiety levels (B) or BLA perfusion 
(D) are shown in warm colors; voxels with negative correlations are shown in cool colors. The scatter plots in C and E are derived from the 
accompanying voxel-wise regression maps shown in B and D and are presented for the purpose of illustrating the range of values only. Data are 
displayed at a threshold of p<.05. The clusters indicated with arrows in B and D met a cluster-wise correction (EWE, p<.05) within the ventral mPEC. 
The peaks of the clusters in B (4, 2, —7) and D (2, 4, —4) were localized to the posterior-most portion of the SGC (with both clusters extending into the 
hypothalamus) using two independent atlases (the Talairach and Tournoux Stereotaxic Atlas [46] and the Wake Eorrest University (WEU) PickAtlas 
[47]; see Methods). Prior work further supports this localization; previously reported sites that have been localized to the SGC (BA25), as well as an 
architectonic mapping of BA25 [68], overlap with the two clusters reported here, with nearby peaks: 4, 2, —4 [69]; -2, 6, —6 [8]; -2, 8, -10 [70]; -3, 
9, -6 [71]; -4, 9, -12 & 2, 11, -7 [72]; 0, 8, -16 [73]. BLA, basolateral amygdala; EC, functional connectivity; Hy, hypothalamus; SGC, subgenual 
cingulate gyrus; mPFC, medial prefrontal cortex. 
doi:1 0.1 371 /journal.pone.0097466.g003 



relationship between BLA perfusion and anxiety levels, was not 
significant. Thus, these data provide further support for the 
Hughlings Jackson, top-down "release" model of negative affect- 
demonstrating associations among all three components of the 
model and providing preliminary evidence for a causal relation- 
ship between poor mPFC-amygdala connectivity and elevated 
amygdala activity and anxiety. 



It should be noted that the mediation analysis in the current 
study was conducted using data extracted from the regression 
analyses; thus the findings only apply to this dataset. It will be 
necessary to repeat these analyses in an independent sample, and 
in clinically ill samples, to confirm these results. Also, the fact that 
the direct effect of BLA-mPFC connectivity on anxiety showed a 
trend towards significance after BLA perfusion was taken into 
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Figure 4. BLA perfusion mediates tKie relationsKiip between BLA-mPFC connectivity and anxiety levels. Mediation analyses revealed 
that BLA perfusion levels mediate tiie relationship between BLA-mPFC connectivity and anxiety (A: iVlodel #1). In contrast, BLA-mPFC connectivity 
did not mediate the relationship between BLA perfusion levels and anxiety (B: Model #2). Values are unstandardized regression coefficients reflecting 
the direct (paths a, b, and c') and total (path c) effects of each relationship in the mediation model. BLA, basolateral amygdala; mPFC, medial 
prefrontal cortex; *p<.05, **p<.01, tp<.001, ttP<-0001. 
doi:1 0.1 371/journal.pone.0097466.g004 



account as a mediator suggests that additional factors, other than 
amygdala activity, influence the relationship between BLA-mPFC 
connectivity and anxiety. Follow-up studies can also consider the 
effects of other potential mediators on these relationships. 

A model of the function of this pathway has been derived from 
findings of studies conducted in rodents. The infralimbic (IL) 
cortex in rodents sends glutamatergic inputs to the GABAergic 
intercalated cells of the amygdala, which, when stimulated, inhibit 
the output neurons of the central nucleus of the amygdala [60,61]. 
Via this pathway, IL is able to inhibit, in a context-dependent 
manner, the fear-related autonomic responses generated by the 
central nucleus of the amygdala. There is some evidence from 
neuroimaging studies in humans that the human homologue of 
rodent IL is the SGC [62,63], although some studies suggest that 
the homologue is located more rostraUy in the mPFC [64] , or that 
there is both a rostral and SGC portion of the mPFC that have a 
role in fear inhibition [62]. 

In the current study, resting BOLD activity of an extended 
region of the ventral mPFC was positively coupled with resting 
BOLD activity of the BLA, consistent with previous findings 



[23,24] and d()S(;ly mirroring the anatomical distribution of the 
amygdala-mPFC pathway seen in non-human primates [25,26]. 
The portion of the BLA-mPFC circuit that was found to be weaker 
here (i.e., showing reduced BLA-mPFC coupling) in subjects with 
higher anxiety or BLA perfusion levels was the BLA-SGC 
pathway. In light of the animal literature and human imaging 
data indicating that fear inhibition is in part mediated by the SGC 
[62-64], our data suggest that the function of projections from the 
SGC to the inhibitory intercalated cells within the amygdala may 
be compromised during negative affect states, leading to elevated 
firing of the amygdala (and elevated perfusion) and increases in 
autonomic arousal and anxiety. 

In the present data, the portion of the amygdala showing the 
greatest elevation in perfusion in relation to anxiety levels 
corresponded roughly (based on an adas-based segmentation of 
the amygdala) to the BLA, rather than the central nucleus of the 
amygdala, suggesting that disruption of the SGC projection to the 
BLA (which is larger than the central nucleus in humans) may lead 
to its cUsinhibition, resulting in over-stimulation of the central 
nucleus and elevated arousal and anxiety. Since the GABAergic 
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intercalated cells reside between the basolateral and central nuclei 
of the amygdala and influence the activity of both nuclei [65], this 
is one of several plausible mechanisms that could account for our 
results. Additional parallel work in humans and rodents will be 
needed to further test this possible interpretation of our fmdings. 

In summary, these results show that ASL is a sensitive tool for 
objectively measuring the neural effects of subde variation in 
anxiety levels. Given that ASL scans are non-invasive, brief and 
have excellent test-retest rehability, ASL could be used in the 
future as a clinical procedure which aids in the assessment of risk 
for future devekipment of psychiatric illnesses, including severe 
anxiety disorders, such as PTSD, depression or other illnesses 
associated with impaired functioning of medial temporal lobe 
structures, such as schizophrenia. Secondly, we also found 
associations among individual variation in amygdala-mPFC 
functional coupling, amygdala perfusk)n, and anxiety levels in 
the same healthy subjects, with mediation results supporting the 
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